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Abstract 

The configurational distribution function, solution of an evolution (diffusion) equation 
of the Fokker-Planck-Smoluchowski type, is (at least part of) the corner stone of polymer 
dynamics: it is the key to calculating the stress tensor components. This can be reckoned from 
[1], where a wealth of calculation details is presented regarding various polymer chain models 
and their ability to accurately predict viscoelastic flows. One of the simplest polymer chain 
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idealization is the Bird and Warner's model of finitely extensible nonlinear elastic (FENE) 
chains. In this work we offer a proof that the steady state configurational distribution equation 
has unique solutions irrespective of the (outer) flow velocity gradients (i.e. for both slow and 
fast flows). 

Keywords: FENE dumbbell chains; Fokker- Planck- Smoluchowski equation; existence and 
uniqueness of solutions; slow and fast viscoelastic flows; Krein-Rutman theorems; 

1 Introduction 

The viscoelastic flow behavior of polymeric hquids is strongly influenced by the complexity of 
various inter and intra molecular interactions. At microscopic level, long chain entanglements are 
a consequence of chain connectivity and backbone uncrossability due to intermolecular repulsive 
exclusive volume forces. Macromolecules diffusion (and conformational relaxation) is slowed down 
due to hydrodynamic drag and Brownian forces. 

Bird, Curtiss, Armstrong and Hassager, together with their collaborators (see [1] and ref- 
erences cited therein), enriched signiflcantly Kirkwood's early ideas [2] and produced a general 
kinetical theoretical framework for both diluted and concentrated polymeric systems. Here, the 
macromolecules are modeled as freely jointed bead-rod or bead-spring chains. One of the simplest 
version of this chain model is the (now popular) Bird - Warner's elastic dumbbell chain, that 
consists of two beads connected by a Finitely Extensible Nonlinear Elastic - aka FENE - spring. 
The salient features of this model, of relevance to this work, are briefly reviewed below, for sake 
of clarity. 

Let X G M", n = 2,3, denote the (microscopic) dumbbell connector vector, y G M" the 
(macroscopic) Eulerian position vector. In the absence of inertia and of external forces, the 
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balance of hydrodynamic, Brownian and intermolecular forces results in the so-called Fokker- 
Planck-Smoluchowski with Dirichlet boundary condition for the configurational function 
ip{t,y,x) diffusion equation. In dimensionless form it looks: 



^ 2De 2De ^ 



over 5(0, 5) (la) 

^\dBiO,5) = (lb) 

In the above equation, -8(0, 6) is the open ball of radius 6 centered at 0, De is the Deborah number 
and = (Vyu)"^ is a tensor which represents the (macroscopic) velocity gradient; the corresponding 
term accounts for the flow type. The second term in the rhs represents the statistically averaged 
Brownian force due to thermal fluctuations in the liquid. The last term, F^'^\ is the elastic force 
that accounts for the dumbbell's elastic response to strain input, for which Warner [3] proposed 
the following expression (valid for ||x|| < 6, with 6 a polymer depending parameter): 

F(^)(i) = ^^-^ (2) 

The above is commonly called the FENE force. Now, as an aside, the model is quite flexible in 
that it may sustain other types of elastic forces: e.g. Peterlin's force (actually a linearized version 
of eq.([2])) usually referred to as FENE-P (see [ll[5]): 

F(^) (i) = ^ ^ = ! ^ (3) 

1- < ||x||2 > 1- <tr(i®i) > 

Asymptotic solutions to the diffusion equation are known for some steady state flows: see (for 
concise presentations see [6], [71 [H [9], [10] ) . They were obtained through series expansions about the 
(known) equilibrium function ipeqi'^)- 

Next, let (t,y) G IR+ x (Q C MJ^). The momentum balance equation reads (see [IT]): 

— + (u ■ Vy) u = uAu - VyP + Vy ■ S, ovcr R+xQ (4a) 
Vy-u = 0, overM+xg (4b) 



where u > 0, u = u(t,y) G M", p = p(t,y) G M. S(t,y) G Sym(M) is the symmetric extra stress 
tensor given by ([!]): 



S(t,y) = /i 



/ F(i)^(t,y,i)di- / ^(t,y,i)dil 

Jb(0.5) Jb(0.5) 



(5) 



UBCO,^) JB{0,5) 

where /x > is a fluid related parameter (actually a given constant). 

One observes that whenever the velocity gradient is such that dui/dyj = aij = constant, 
'^iCii = 0, and is a solution of (fTal) - (fTb|) . then S defined in equation ([5]) is always independent 
of y, hence Vy -8 = 0. In such a situation there exist u and p so that (Hall -( ]4bl) are solved. That 



this is indeed the case may be inferred from the following. Using Einstein's summation convention 

du ■ 

over dummy indices, Ui = aijyj + c, therefore ^^fc = (^ik [^kjUj + c] = aikakjVj + di, di = c^-,^aik. 
Hence [Vy • (u ® u)],. = aikakjUj + di = Cdji/j + di. Therefore Vy ■ (u eg) u) may be expressed as 
Vy ■ (u ® u) = — VyP, where p = — (1/2) aijUii/j — dii/i, since the matrix of entries a^j is symmetric. 
We conclude that for any traceless matrix A whose entries aij are constants, and for a steady 
state, homogeneous flow solution ip{x.) - i.e. independent oft and y - to equations f|Taj) - ffTbl) . there 
exists a steady state solution to fHa|) - fj4bl) given by: 

Mi(x) = aijUj + c (6a) 



^aikakjUiVj - c ^ aikVi (6b) 



P = 

i,k 

and with S given by eq(l5]). 

For this work we shall consider u as being given by eq fl6al) . where A is a given matrix, and we 
shall prove the existence of a solution t/^, independent of t and y, to f|Ta|) - ffTbl) . 

Before proceeding further, we pause for the following important observation. The solution 
ip to (fTaI) - (fTb|) we inquire about - being a probability density - has to be non-trivial {ip ^ 0), 
non- negative and integrable. As ^/^ = is a solution to the aforementioned problem and as we 



have to mind about non-trivial ones, the solution non- uniqueness must be compulsory. Next, we 
know from fi2[ [13] that ip = is the unique solution to f|Taj) - fllb| ) whenever F'^'^^ is an element of 
L''"{B{0,6)), r > n. Therefore, what makes possible the existence of non-trivial solutions, is the 
fact that F('=) is NOT an element of U{B{{], 5)), r > n (in fact F^'^) is not an element of ^{3(0, 6)) 
for any r > 1). 

Now, as FW = Vxt/(i), with U{±) = - (5^/2^ log (l - ||i||V5^), equation 1^ is usually 
re- written as (see [H]): 



2De 



M(i)V, 



M 







where the function M : 5(0, 5) ^ M is given by: 



1 / UP^'/' 
M(5c) = i(l-& 



where J is a normalization constant so that: 



(7) 



(8) 



M (±) di = 1. 



B{0,5) 



(9) 



Next, for sake of generality, we replace 0x by an arbitrary function k : B{0,6) M", and the 
problem we focus on can be re-formulated as: 



1 



2De 



V- 



MV 



M 



^1 



dB{0,5) 



where a 



ip (x) dx = a 

IB{0,5) 

is a given constant. 



(10a) 
(10b) 
(10c) 
(lOd) 



meas((5) 

Next, for notation convenience, we carry out the variable change x = x/5. This transforms the 



domain n into n = 5(0, 1) = {x G M", ||x|| < 1}. Let us denote ^(x) = ^(x), k(x) = 25Dek(x), 



6 = P/2 and let 



M : Q 



i, M(x) = (1 - ||x 



\2\5 



Then, equations fllOap - (]10d|) become in Q: 



-V- 

^\dn = 



M 



V ■ (k^) = 0, Vx G ^] 



V^(x)dx = b 



(11a) 
(lib) 
(11c) 
(lid) 



with b > and k : i7 — M" given. As in practical situations S is (roughly speaking) close to 10, 
then 6 is close to 50. 

The goal of this paper is to prove the existence and uniqueness of a solution to the system of 
equations flllal) -( fTTdl) . We easily see that the aforementioned problem can be also formulated as 
following: prove that is a simple eigenvalue of the operator (denoted from now on L) defined 
by the Ihs of (II lap and the boundary condition (lllbp . with a corresponding non-negative and 
integrable eigenvector. In fact, we will prove that is the principal eigenvalue of L in the sense 
that the real part of any other eigenvalue of it is non-negative. To achieve this we use the classical 
Krein-Rutman theorems, in both weak and strong senses, on an appropriate operator obtained 
from L. This will also entail that ip is positive over Q and behaves like M on the boundary dQ. 

The boundary value problem problem with unknowns u and ^ as presented in (diD, ([TbD, dUD, 



(Hbj), (|5]) has attracted the attention of several investigators working in the area. For instance, in 



[15] Zhang and Zhang proved the existence of a local in time, regular solution to the system formed 
by equations ffTal) . ffTbl) . (Hal) , fjlbj) . and (jS]). The existence of a global in time solution was proved 
in [16] by Lin et al, and that in a particular case referred to as the "co-rotational" velocity field, 
that is, in equation (fTal) the term 6 = (Vm)"^ is replaced hj 6 = VyU — {VyuY- Moreover, for this 
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same system of equations, in [T7] , Barret et al offered a proof for tlie existence and uniqueness of a 
solution to a regularized problem associated to the system fITal) . fflbl) . fj4a|) . fjlbj) . ([5]). Next, in [T8] 
Lelievre et al proved the existence and uniqueness of a local in time solution to the one dimensional 
motion system of equations in which the Fokker-Planck-Smoluchowski equation is replaced by a 
stochastic diffusion differential equation. 

In [19] Du et al focused on the Fokker-Plank-Smoluchowski evolution equation only, assuming 
a steady and homogeneous macroscopic velocity field. For this they proved the global in time exis- 
tence and uniqueness of a solution. For the corresponding steady state problem, the forementioned 
authors proved the existence of a solution only in the particular case where the tensor 6 in ( ITal) is 
either symmetric or antisymmetric. Some of their numerical simulations suggest the existence of 
steady-state solutions for arbitrary 6. 

In this work we do prove the existence and uniqueness of steady state solutions for arbitrary 

e. 

As an aside, in [T?] Degond et al provided arguments in support of the validity of an asymptotic 
expansion solution, valid for small De numbers, first obtained in [1]. 
This paper is organized as follows: 

• in Section [2] we state the main steady state existence and uniqueness result, 

• Section [3] addresses some important functional analysis preliminaries, 

• Section mis devoted to proving the conclusive existence and uniqueness result. 
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2 Functional framework. Presentation of the main result. 



Let the following spaces be defined as: 



X < oo 



u 

— + M 
M 



M 



dx < oo 



endowed with the norms 



2 \ 1/2 

— dx 



and respectively 



u 

— + M 
M 



M 



dx 



1/2 



(12) 
(13) 



It is clear that L\.j is a Hilbert space. To see that Hlj is also a Hilbert space, let us remark 



that 



Hli with = |t; G Li^iSl), J {Mv^ + M |V^;|^)dx < ooj . 

It is well-known, as being a classical weighted Sobolev space, that Hj^ is complete (see for 
exemple Theorem 3.2.2. (a) in Triebel's monograph [20]) when endowed with the norm 



1/2 



Jn / 
Since the application -ip G Hl^j Mip G Hlj is an isometry, we deduce that Hl,^ is complete. 



For any ip G Hl,j{Q) we denote \(p\i the semi- norm on i/]^ defined by \(p\j := M 



M 



dx. 



Moreover, (if|^(n))' denotes the corresponding dual space and one has the canonical embedding 



We now endeavor to search for solutions to equations flllal) - fllldl) that are elements of Hlj, as the 
trace on dQ of any u G Hl,j{Q) is zero (see also Proposition 13.31) . 

To achieve this goal, equation (11 lap is first multiplied by (f/M, with (f G V{Q) and next 
integrated over fl. It gives: 



MV 



M 



Definition 2.1. ip G H\,j{VL) is a weak solution of the system f ilial) - fllldp . provided that: 



and moreover, flllcp and ( llldl) are satisfied. 



Next, let the operator L : H\^{Vl) {Hlj{Q))' be defined as: 



{L(u 



Now, L is well-defined, due to: 



m) \mJ 7o Vm 



< 



< 



U 

M 



M 



dx+ / Ikl 



U 



Ml/2 



M 



dx 



< Mhi,\Mh\, + M\l^ Mli,\Mhi, 

It is now clear that our problem is tantamount to finding an element G H\j^ such that 



= 



> 



/ '?/^(x)dx = 6, 
Jn 



that is, must by a non-negative and integrable eigenvector of L corresponding to the eig 
0. 



For any P >0, let: 



Xf3:= {ipeC{n) ,3c>0 s.t. |(^(x)| < cM^(x), Vx G 



Xg is a Banach space endowed with the norm 



Mx,:=svip^^^ = mf{c>0 s.t. |<^(x)| < cAf^(x), Vx G ^]} . (19) 
Remark 2.1. Xq = (fi), endowed with the usual norm Hv^llx,, = sup \ ip\. 
Remark 2.2. For Pi < P2, the continuous inclusion Xp^ C Xp^ holds true. 

cont 

Let the cone Pp C be defined as: 

Pp := e Xf,, (^(x) > 0, Vx G fi} . (20) 
It is clear that P/3 is a reproducible cone for the space Xp, that is X/3 = Pp — Pp. 

o 

It can be easily seen the interior Pp of Pp is given by: 

Pp = i^^eXp s.t. mf >o| = {v9GX^, 3ci >0 s.t. (/p(x) > CiM^(x), Vx G 1^} 

= {y? G "^(r^), 3ci, C2; < ci < C2, s.t. CiM'^(x) < (/?(x) < C2M'^(x), Vx G fi}. (21) 

We now state the cornerstone result of this paper: 

Theorem 2.1. [Existence and uniqueness theorem] Let b > 0, 6 >8 and k G (Vr^'°°(i7))" . 
Then there exists an unique solution ip to the system ( fll7ap -( TT7c|) ). Moreover, this solution belongs 

o 

to Pi which amounts to say that ip is continuous in Q, and there exist ci, C2 with < ci < C2 
such that 

CiM(x) < V^(x) < C2M(x), V X G 

Remark 2.3. We assume throughout this paper that 6 and k comply with the hypotheses of The- 
orem \2.1\ Given the physical model under consideration (for which we gave a suitable description 
in the Introduction section), such an assumption does not lower down the level of generality. 
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3 Several preliminary results 



3.1 Basic facts 

In the following we denote for any real a the operator : H]^{Q) — > {H\j{Q))' given by 
L + ald-i Id being the identity operator. 

It is assumed, throughout this paper, that a is large enough so that: 



where 



a > max -||k||i,oo(n) + 1, 4Ao + XqU + 2Ao||k||ioo(Q) + ||V ■ k||ioo(s^) 



Ao = 2(||k||ioc(n) + 1) 



(22) 



(23) 



Proposition 3.1. The operator is invertible. 



Proof. Let / G (iJ^(f2))', arbitrary. We have to prove the existence of a unique solution u e Hj^ 



to 



aa{u,ip) = {f,(p), V<^ e Hi 



M 



(24) 



where, in the above. 



a / — dx 

f2 M 



(25) 



Next, to use the Lax-Milgram theorem, one only needs to prove is coercive as all other theorem 
constitutive assumptions are obviously fulfilled. 

The fact that is coercive is an immediat consequence of the inequality 



];^k,.v(^)dxi</jki 

< 2 + llklli-lbllii^) 



MV2 



Ml/2 V / 



dx 



and of the choice of a. 



(26) 
□ 
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Let then Ba : {Hm)' —>■ Hm denote the inverse operator of La- 
Clearly 

BaE C {(Hlf)', Hm) . 

and, also, 

Ba E C {L\j^ L\j) . 

Lemma 3.1. [Weak Maximum Principle] Let f G {H\,^)' , / > 0, and u = Baf ■ Then u > 0. 

Proof. The proof is classical and consists to choose ip = u~ m. the corresponding variational 
formulation, (see for exemple [21] for the non-degenerate case M = 1). □ 

Lemma 3.2. [Comparison Principle] Let Vt' he an open set such that Vt' C Vl, and let VL' denote 
its closure. Let u,v E Hlj{Q) so that LaU and LaV he functions well defined on Vl' . Assume that 
LaU > LaV on Q' . Then: 

• Case 1. If Q' G Q and if u > v on dQ' , then u>v on Q' . 



• Case 2. If^l — ^l'C^l and if u > v on d{^l — ^l'), then u > v on Vt' . 

Proof. Let w = u — v E Hl,j{Q) and / = La{u — v). Then, for Vy? G Hl^{fl), such that ipln-w = 0, 
one has: 

/ MV (^) ■ V (^) dx' - / . V (^) dx' + a / ^dx' = / ^dx' (27) 
We now take in [27] (/j : J7 M defined by 

{w~ on Vl' 
(28) 
on VL-VL' 

and we easily obtain the result. □ 
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We now introduce for any /5 > 



(29) 



Actually, L2,/3(fi) is a Hilbert space endowed with the norm 



w\\L2,Bm 



(30) 



We have, as a straightforward consequence of Theorem 6.2.5 of [H], the following continuous 



inclusion: 



(31) 



Next: 



Proposition 3.2. (a) If (3 > 1/2 - 1/(25), then Xp C L^. 

cont 



(h) Ifl3> 1/2 - 3/(25), then Xp C {Hi,)'. 

cont 

Proof. (a) Let ip G Xp, arbitrarily. Then, |v5(x)| < M^(x) from which we get that: 



dx < llv^lli 



M 



dx = II (/9 1 



X, 



I (l-||xf) 



2\{2,0-l)<5 



dx 



(32) 



However, (1 - HxlH^^'^"^^^ dx < +oo iff /5 > 1/2 - 1/(25). 



(b) Let if G Xp and G H\j arbitrarily. We have 



^dx 



< \Mx, I M'^-M^ldx 



(33) 



However: 



M'3-i|^|dx< ||M-V2-i/^^||^^^^^ ||^"-^/'"''/1L.(n) 



(34) 



^2,1/2+l/i 



Moreover, the L^-norm of M^-^/^+V-^ is finite iff /? > 1/2 - 3/(25). 
Using also the continuous inclusion (!3T|) we have the result stated. 



□ 
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Proposition 3.3. Let (3 he such that < /3 < 1/2 - 1/5. Then, for any u € Hlj, G HI{Vl); 



moreover, 



u 



u 



< C \\U\\ rrl . 



u 



_^i/2-/3_ actually need to prove that -—pr e H\n). From this, 



Proof. Let v ■■ , • — — — - j^....^ , ,^ 

since M^/^"^ G ^^(H) and M^/^-/?!^^ ^ g, it will follow that v G HI{VL). 
To begin with, notice first that j^pj^ ^ -^^(^)) as u G L^^(n). Next, 



u 



u 



u 



V ( M^/^— ] = M^V ( — ) + V 



(35) 



Now, M^/V G as u G ifl^. Let us next show that V (M^/^) ^ G L^(r]). One 



has: 



u 



VM 



M 2M3/2 MV2+1/5 2Mj-i^ 

eL°°(n) 



(36) 



Next, by ^ 



u 



j\/fl/2+l/<5 



Hi 

G Z/^(ri), from which we infer that 'j^pj^ ^ H^(n), and further on that 



u 



G ifQ(r2). It is easily deduced that: 



u 



MP 



< C ||n|| rrl 



(37) 



□ 



Remark 3.1. Taking f3 = in Proposition \3.3\ we deduce that u G Hq{VL) whenever u G H\j{VL), 
which triggers that the trace of u on the boundary dQ is equal to zero. 

Proposition 3.4. Let if G Xp, /? > 1/2 + 1/(25), be such that Vy? G (X^)" with 
7 > 1/2 - 1/(25). Then ^ e Hj^. 



Proof. From Proposition 13.21 a) we have ip G L\j. Next, 
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/ M 




2 

dx = 


/ M 


Jq 






Jn 



VM 



dx 



However, 



(38) 



[ — |V(/?|Mx<c /" — M2Tdx<+oo 

provided that 5(27 — 1) > —1, which amounts to7>l/2 — 1/ (25). 
Next, 



(39) 



M3 



-dx < ci 



M3 



-dx < C2 



< cJ {1 
Jn 



(I -Mi" 



(1 — ||x|| j dx 



I ,,2x2(5/3+2(5-2-35 , 

|x|| j dx 



(40) 



For the above equation fHOj) to hold true it is necessary that 25(3 — 5 — 2 > —1, i.e. (3 > 
1/2 + 1/(25). 



□ 



Since G X^g for any /? > we have the following direct consequence of the above result: 
Proposition 3.5. For any /3 > 1/2 + 1/(25), we have M^^ e Hl^. 



3.2 Continuity and compactness properties of 

The goal now is to appropriately introduce several {¥', Y") paires of Banach spaces such that 
is well defined and continuous from Y' to Y". Some compactness properties of Ba, needed further 
on, are also proved. 
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Lemma 3.3. (i) Let (32 e R such that 1/2 - 3/(25) < (32 < 1/2 - 1/6. Then G C{Ll^, Xp^) . 
Moreover, is a compact application from L\.j onto X^^ ■ 

(ii) Let (3u(32 e R such that 1/2 - 3/(25) < (32 < 1/2 - 1/5, and (3i > (32- Then B^ e 
C{Xi3j^, Xp^) . Moreover, B^ is a compact application from X^^ onto Xp^. 

Proof. The proofs for the above two statements are pretty much similar in nature; henceforth, we 
offer below a global proof, and pause wherever necessary to particularize it so to get the results 
in either (i) or (ii). Keeping that in mind, let / G L\j (for the (i) part) and / G Xp^ (for the (ii) 
part) and let u = B^f . Observe that: 

Ll C {H\,y (41) 

cont. 

and that 



(42) 



as consequences of Proposition I3.2[ In both cases / G {H\i)' and u G H\j solves the equation 



- V- 



MV 



u 
M 



V • (ku) + au = f, u e Hi 



M 



(43) 



We also have 



\u\\hI < c||/||l2 for the part (i) 



and 



\Hlj ^ ^\\J lix^j for the part (ii). 



ImII TTi < c 



Denote f (x) = u (x)/M^2(^x); we first take on to prove that v is bounded on fi, which prompts 
that u belongs to X^^. 

Making use of the fact that u = wM^^ j^to leads to: 
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-V • [M^'iVv) + {p2 - l)M^'-\VM)v] + (V ■ k)M^^t; 
+p2M^^-\VM) ■ kv + kM^2 ■ (Vv) + aM^^v = /, Vx G 

which, after a few re-arrangements, can be re-written as: 



Av = g 



where 



/ 



+ 



(2A-i)^-k 



M 



Vv 



+ 



AM 



(A-l)=^ + (ft-l)^^-V.k-,8.^.k-a 



We also deduce from Proposition 13.31 that 



In the following we shall obtain some convenient estimates for the function g. We have 



Vv = V(— 



rl-/32^ 



- M^-^^V (^^) + V(M 



GL2(f7)as«em, 



Using equation (l46ll we get: 



where, in the above. 



91 W 



52(X) = (1-/32 
1 



M 

VM 

(2/52-l)^-k 



VM 

M^2+i 



+ 



M/32 



AM IVMP VM 

(/?2-l)-rr + (/32-l)^r7^-V-k-/?2^-k-a 



M 



M2 



M 
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For the (i) part of Lemma 13.31 one has: 



/ 



L2 



(50) 



while for the (ii) part of Lemma 13.31 one gets: 



/ 



Moreover, 



VM 



L2 



L°° II'' 11^, 



/3l 



(51) 



1 



1 



M llxHi 1 - ||x||2 MV-s" 
Therefore, the above leads to gi G L°°(fi). We then deduce giM^/'^S/{u/ M) G L'^{Vl) and 



\giM^/^V{u/M)\ 



(52) 



Now, observe that: 



which imphes 



5^2 (x)m 



lIxlHl MV2+1/5 



^1/2-/32-1/5 

V ' 

eL°°,as/32<l/2-15 



We deduce with the help of inclusion f l3T|) that 



|5'2M||L2(n) < ci||m||hi 



(53) 



and further on, from (|49l), ([50]), (Ell), ([52]) and ([53]), that 



(54) 
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for part (i), and 

hh^n) < c||/|U,, (55) 

for (ii) part. 

Now since v satisfies ( H6l) and ( H71) we obtain v G H'^{fl), and ||f ||_f/2{Q) < c||5'||L2(n). By the 
Sobolev's inclusion C C(i7), n = 2, 3 it follows that f G C(i7), and that Utill^/f^-, < ci||5f||j^2(f^). 

cont ^ ' 

Thus u G X/32 and 

\\u\\x,^<cMLHn) (56) 
Next, making use of (IMI) and (15^ . we deduce, for part (i) that G =Sf (L|^, X^^), and for part 
(ii) that Bo, G ^(Xfs^Xp,). 

In order to show the compactness of Ba, let {fg)gem be a bounded sequence in L\,j for part (i), 
and in X^^ for part (ii), respectively. Denote Ug = -B«(/g) G X^gg, and Vg = Ug/M^'^. Next it is 
proved that Vg is bounded in H'^iVt). As the domain VL is bounded, the inclusion H'^iVt) C C(f2) 

is compact; hence there exists a subsequence q' of g and a f G C(f2) such that Wg/ -j^ v. Denoting 

c{n) 



u = vM^^ , we have that u G Xp^ and sup — ^^^^^ ^ ^ • Therefore m^' u with respect 

5'— ►+00 



to the X^2 space topology. 



□ 



For any r > 0, let us denote Vtr := {x : ||x|| < r} = 5(0, r). 

Lemma 3.4. Let /3 be such that 1/2 + 1/(25) < (3 < 1. Then B^ G ^(X^_2/5, X^). 

Proof. Proposition O b)gives that X^_2/5 e (ifj^)' (as 1/2 - 3/(25) < /3 - 2/5), which entails 
that the operator B^ is well defined over Xp_2/s- Let / G X^_2/5 and u = Ba{f). We have to 
prove the validity of the following assertion: 
There exists A' > independent on / such that 

AM^i^) - m(x) > 0, Vx G 1) (57) 
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and 

y4M^(x) + u(x) > 0, Vx G Q (58) 

where we denoted A = A'||/||/3-2/<5- We shall provide a proof for the first one only, i.e. for (!57|) . 
as the other may be proved similarly. The proof for (1571) relies on the Comparison Principle stated 
in Lemma [3.21 



One has: 



L„(M^) = M^(a + V-k) + M^~^[/3k-VM+(l-/3)AM] 

- M^-2(l -/5)2|VM|2 (59) 



As: 



VM = -25x(l - ||xf )'^"^ 

AM = -2n5(l- ||xf)^-i+ 4(5(5 -l)||x||2(l- ||x|n^-2 

then, 

L^M^) = ao(x)(l - ||x|n^/5-2 + ai(x)(l - M'Y^-' + a2(x)(l - Hxi^^^ where 
ao(x) = 45(l-/5)(5/5-l)||x|p 
ai(x) = -[25/3x-k + 2(l -/5)n5] 
a2(x) = a + V ■ /c(x) 

It is clear that ao(x) > 0, V x G fi. Next, since / G X^_2/5, we deduce 

-/(x) > -||/IU,_,/,M/^-2/^(x) = -||/||x,_,/,(i - M'Y^-' 

Then: 

AL.(M0 - / > [Aao(x) - ||/||x,.,/J(l - ||x|n^^-^ + Aai(x)(l - M^f^-' 
+ Aa2(x)(l-||x|n^/^ 
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In the following, we restrict ourselves to ^2 — r2i/2, henceforth ||x|| > 1/2. Then 



ao(x) > a° := 6{1 - (3){6(3 - 1) > 0, Vx G - 



and 



AL,(M^) - / >Aa°(l - ||xf 



+ ^(i-iixr) 



a, 



(60) 



02 (Xj 



^(1- ||X||^)^], VXGI]- 1^1/2 



Assume f ^ (this is not too restrictive as, whenever / = 0, the inequality (1571) is satisfied with 
A = 0). 

Let us choose tq G [|, l[ close enough to 1 such that 



1 — 1 

^sup|ai(x)| < - 

^ o^sup a2X <- 



(61) 



On the other hand, from Lemma 13.31 (ii) with f3i = f3 — 2/6 and /?2 = min{/5 — 2/5, 1/2 — 1/6} we 
deduce m G Xjj^ (since / G X^J and 



u\\c{n) < Ci||m||x^^ < C2||/||x 



(62) 



Take now 



A = max 



C2 



X 



13-2/S- 



Clearly, from (160D. (16T]) and (162|). 



L«(AM^ - u) = AL„(M^) - / > on - a 



and 



AM^ > u on Q 
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respectively. Invoking the Comparison Principle (Lemma I3.2p and the fact that u, G leads 
to 

AM^ >u on Q-Qro, which implies AM^ >u on VL. 
This ends the proof. 

□ 

Lemma 3.5. G ^(Xi_i/5,Xi) 

Proof. Let / G Xi_i/5 and u = Ba{f). Let : ^ M, W{y.) = e^H'^H'. The job is now to prove 
that there exists A > and A' > independent on /, such that 

|u(x)| < AW^(x)M(x), VxGl] 

where we denoted A = y4'||/||xi_i/^-- Actually we take on to prove AWM > u only, as AWM > —u 
can be proved similarly. To achieve this, we again make use of the Comparison Principle. 
One has: 

L^{MW) = M[-AW + k ■ VW + (a + V ■ k)W] + VM ■ {kW - VW) 
However, VW = 2\xW, AW = {2Xn + \'^\\xf)W, hence: 

Lc,{MW) = Miy(-4A2||x||2-An + 2k-Ax + a + V-k) 

+ 25Mi-i/V(2A||xf -k-x) (63) 

Let us take A = Aq with Aq given in (l23l) . We obtain 

2Ai|x||2 - k-x > 1. 

From hypothesis (122|) on a we obtain 

- 4A^||xp - An + 2Ak ■x + a + V- k>0 on Q - Q1/2 
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which gives 

La{MW) > 26WM^~^^\ Vx G 1] - l^i/z- (64) 
Next, as / G Xi^i/s, one gets 

-/(x)>-l|/l|x,_,/,M^-V^(x),VxGfi 

and invoking further on (IMll leads to: 



AL^{MW) -f> [2SWA - \\f\\x,_,jM'''^', Vx G 1] - VA > 0. (65) 
Choose A > such that 26AW{'x) > Vx G ^2 — f2i/2- For instance, any A such that: 



^>^II/IU,_,/, (66) 

will fit in. Then: 



La{AMW-u)>0 on Q-Qi/i (67) 



On the other hand, one needs to choose A so that AMW > u holds true over ^1/2- We proceed 
as in the proof of Lemma 13.41 A sound choice for A is one such that 



A min W{^) min M(x) > maxM(x) (68) 

Next, min W{^) = 1, min M(x) = (3/4)'', and we are left over to inquire about max'u(x). 
To get an answer to, we shall call in Lemma [3.31 with /3i = 1 — 1/5 and /32 = l/2 — 1/5. One has 

\\u\\c(n) < ci|k||xi/2_i/, < C2||/|Ui_i/,- 



23 



Then, one may choose A > (4/3)'^C2||/||xi_i/5 t° ensure (l68ll holds true. Finally, taking into account 
(l66l) . we are left to choose 

A = max{l/(25),(4/3)%}||/|U,_,/, 

and we end the proof exactly as in Lemma [331 taking into account the fact that MW G if^/, so 
that the Comparison Principle can be made use of. 

□ 

3.3 Strong Maximum Principle for the Ba operator 

This section aim is to prove the following "Strong Maximum Principle " property for B^. for any 
/ G Pi - {0}, BJ eh 

The following weaker result is first proved. 
Lemma 3.6. Let f e Pi, f and u = B^f . Then ■u(x) > 0, Vx G Vt. 

Proof. We adapt here the classical proof for the case where M is equal to 1 (the non-degenerate 
case; see for example Gilbarg and Trudinger [K] or Evans [21]). 

We remark first that u is continuous on Vt. Assume 3x G i7 such that m(x) = 0. 

Denote Ve := {x G m(x) = 0}, K := {x G VL, m(x) > 0}, V; U K = ^■ 

By hypothesis VJ= 7^ 0, as well as K 7^ 0. It is clear that Vg is open and that dVe ^ dVL. Let 
zo G 9Ve n 7^ 0; then m(zo) = 0. Denote d = inf |zo - z| > 0, thus |zo| = 1 - rf. Let ri g]0, rf/4[ 
be small enough, and fix xq G Vg such that |xo — zo| < ri. As is an open subset, there exists 
r2 > such that i?(xo,r2) C Vg. Therefore choose Vq = sup{r s.t. i?(xo,r) C Vg}. Then there 
exists yo G i?(xo,ro) fl K 7^ 0, u(yo) = 0, and |yo — Xq] = Tq. This prompts Tq < ri, hence 
one may choose a small enough tq. Thus, u{yo) = 0, m(x) > 0, Vx G -B(xo,ro). Let the function 
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w be such that w : S(xo,ro) — > R, u'(x) = e '^H'' ^""^ — e '^^o, where A > will be later chosen 
conveniently. Denote also by w the continuous extention of w at on Q and note that w e H\j{Q), 

w(x)|xGai3(xo,ro) = 0, and w(x)|xeB(xo,ro) > 0. 

Next, we take on to prove that 3^4 > such that w(x) > Aw{yi), Vx e -B(xo, ro) — -B(xo, ro/2). 
To achieve this we shall make use of the Comparison Principle. We actually evaluate L^iu — Aw) — 
f — ALa{w), / > 0, and prove that La{w) < for any Vx e -B(xo, ro) — -B(xo, ro/2). 



Basic calculations lead to: 



(69) 



where, in the above 
E2 



-Aw + — — + k -Vw 
M 



V- — — + V-k + a 

M 



w 



AM |VMp 



M 



M2 



+ V • k + a 



Using the expression of w we find 



El 



4A^|x-xoP + 2An + 



45 Ax ■ (x — xq) 



- 2Ak • (X - Xq) 



-A|x-xo|2 



and 



Eo — 



2nS 45((5-l)||x| 



1- ||X||2 

+ V • k + a 



+ 



4(52||x||2 



(1-||X||2)2 (l-||x||2)2 



(70) 



(71) 
(72) 



Next, observe that ||xo|| < || < 1 — (i+ri, so for any x G i?(xo, ro), one has ||x|| < 1— ci+ri+ro < 
1 - d/2. Therefore 1 - ||x||2 > d{l - d/A) > 0. Denote d^ = l/(d - d^/A) > 0; hence 

1 



< do, VxeS(xo,ro) 



(73) 



25 



Then 



El < [-AVq + 2\n + 45Arorfo + 2A||k||Looro] e^^l''"''"!', Vx such that ^ < ||x|| < tq. 



We also have 



E2 < \E2\ < [2n5dQ + A8{5 - l)dl + AS'^dl + || V • k||i=o + a] e"^!^-^"! 



which imphes 

Laiw) < + A(2n + 46rodo + 2||k||iooro) + 2n6do + (85^ - A6)dl + || V ■ k\\L^ + a] e"^!™!', 

Vx such that < ||x|| < tq. 

(74) 

Then one may choose a A > large enough (with A depending on zi^ and ro) so that the rhs be 
negative, i.e. LaiyS) < 0. Therefore 



L^iu-Aw)>0, Vxe5(xo,ro)-5(xo,ro/2), > 0. 



(75) 



Next, as m(x) > in i?(xo,ro) and m is continuous in Q , one has inf ^^(x) > 0. Choose 

xeaB(xo,ro/2) 

^ > such that inf m(x) > A sup w(x) = A ( e"^"""/"^ - e~^''o ) . Then choose 

xe9B(xo,ro/2) x6aB(xo,ro/2) ^ ^ 



A 



inf 'u(x) 

xeaB(xo,ro/2) 



/ 



, to get 



u > Aw on 9i?(xo, ro/2). 



We also have 



u > Aw on (9i?(xo,ro) 



Then the inequality (175|) and the Comparison Principle give 



M > Aw, Vx G 5(xo, ro) - 5(xo, ro/2). 



(76) 
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Next, the interior regularity property gives u G C^iVl). Let v = — denote the outward normal 

To 

vector at yo G 5(xo, ro). Then ^(yo) = lim -u [yo - s(yo - xq)] . 

ou rn s-^o S 

With the help of inequality ( (!76|) ) it easily follows that: 

l^(yo) < -2\Arle-'^o < q. (77) 

On the other hand now, yo is an interior point at which u reaches a minimum ('u(x) > on f2, 

Ou 

''^(yo) = 0); this entails Vu(yo) = 0, hence ^(yo) = 0, which contradicts inequality fl77|) . This 

last argument ends the proof. □ 

The main result of this section is 
Lemma 3.7. 

5„(Pi - {0}) cA . 

Proof. Since Lemma [3751 gives Ba{Xi) C Xi, it suffices to prove that for any / G Pi, / ^ 0, there 
exists c > 0, such that 

m(x) > cM(x), Vx G (78) 

where in the above u = Ba{f). 

The difficulty here is to lower bound m in a neighborhood of the boundary of Q. 

Let Wo : ^ Wo(x) = e-^"^"' - e"^, with A > 0, and 
W^i : ^ M, W^i(x) = [Wo{^)]\ Observe that Wo(x) = e'^ p^^-"''"') - ij ; using now the 
inequalities 

-1> z and e"" - 1 < ze^ with z = X{1 - ) 

we deduce 

Ae-^(1- ||x|n < |W^o(x)| < A(l- ||xf), VxGfi (79) 
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which imphes 

A^e-^^M(x) < |iyi(x)| < A''Af(x), Vx G fi. 
Then Wi e Xi which triggers Wi e L\[. 

The foUowings hold true as weU: VW^o = -2Axe"^ll''ll' and VWi = -2X5xe'^^^^^^^WQ'\ Inequahty 
(1791) leads to |Viyi| < cM^'^/^. Finally, Proposition 1331 gives Wi G Hl^. 

We take on to proving La{u - AWi) = f - ALaiWi) > on i7 - fi^, where A> and r] g]0, 1[ 
will be chosen later. As / > 0, we need to prove that La{W) < 0. One has: 

Carrying out the calculations by making explicit VM, AM, etc, leads to 



L^{Wi) = 2A5ne-^IWlVo^-^ - 4A25||xf e-^""" Vq^^' - AX^6{6 - l)||x|| V^H^HVo^-' 



+ 4X6 \ ""^ „, e-"IWlVo^-^ - 2X6k ■ xe-^H'^H V.'-^ - ^ , (80) 



II 2 II II 2 

^" "-S /ir2 ll^ll w5 , ^V7 1, , „AW<5 



+ ^^^^ " ^ V-||x||T ^° " T^RF^° + (V ■ k + a W 
Denote now y = 1 — ||x|p. Expanding about "y close to 0" leads to 

e^i^ with zi G [0, A]. 



Wo = Xe-^y il + ^y ' ^^y^ 



= A^e- V (^1 + ^ + y'Hy, A, 7) J (81) 

where, due to the fact that ?/ G [0, 1], h is such that 

\h{y,X,'y)\<h 

where h is a positive constant depending on A and 7. Next, for any 7 > 0, 
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2 6 
We then have, for any 7 > 



2) 



with Z2 E [0, A7]. 



Expand the right-hand side of ( !80l) in power series w.r.t. y. Using ( ISTl) and ( l82l) and taking 
into account the equahty = 1 — y, one remarks that the coefficients of the leading term y^~^ 
vanish, so after some lengthy (and awkward) algebra one gets: 

LaiWi) = ai(x,A)/-i + a2(x,y,Ay (83) 
with (84) 
ai(x,A) = 2AV'^^ (-^^A + - (5x • k) (85) 

and 02 a function satisfying 

|a2(x,?/, A)| < 02 
where 0,2 is a positive constant depending in A. Next, 

— (5^A + 4(5^ — (5x ■ k < — (5^A + 45^ + (5||k||ioo and with a suitable choice for A, such as: 

A = ^(452 + <5||k||i^ + l) (86) 
one gets 46^ — 5^A — 5x-k< —1, which gives 

ai(x,A) < -2AV^\ (87) 

Therefore 

UW,) < -2X'e-'V-' + a2y' = -2X'e''V-' (l - |A"^e^^y) . 

Then one may take y small enough (i.e. ||x|| close to 1) such that La{Wi) < 0. 
It has thus been proved that 3?7 g]0, 1[, close to 1, such that 

Laiu-AWi)>0, VxGl]-fi^, WA>0. (88) 
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Next, from Lemma [3.61 we have m > over Q. Since u is also continuous, A may be chosen such 
that min m(x) > AmaxWix). Take A = min u(x). Such a choice leads to 

< 2* 

m(x) > AWi{x), Vx e (89) 

Since u, Wi G Hj,^, use of equations (1551) . (15^ and of Comparison Principle allows one to infer 
that m(x) > yliyi(x), Vx G r2 — Qri- One more use of fl89l) implies that, in fact, this inequality 
holds true on the entire fl domain. Now the inequality (1751) gives the result. 

□ 



4 Proof of the main result 

We are now in position to give the proof of Theorem 12. 1[ 
Step 1. 

From part {ii) in Lemma 13.31 one infers 

Ba G Xi/2-l/<5)- 

From Lemma [3.41 we obtain 

Ba G =^(Xi/2+(2j-i)/5, ^i/2+(2j+i)/5), for any j G {0, 1, ■ ■ -jo} 

where 

r 127 + 1 

Jo = max j J G N, - + < 1 

(in other words, jo is the unique natural number belonging to the interval [| — |, | — |[)- 

Due to the inequality | + ^^^y^ > 1 — | we have the inclusion Xi/2+(2io+i)/<5 ^ ^i-s/s- Using again 

Lemma [3.41 we obtain 

G ^(Xi, Xi_i/,). 
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Finally, from Lemma [3.51 we deduce 

where the following notation has been used: 

Sc. = Bi'+\ 

As Ba is a compact operator that maps Xi onto Xi/2-i/s (see part {ii) in Lemma 1373]) it follows 
that Sa is compact as well. 
On the other hand, Lemma 13.71 gives 

Sa{Pi - {0}) CPl . (90) 

We are now in a position that allows to make use of the strong version of the Krein-Rutman 
theorem (see for example [22]) to the operator Sa, the Banach space Xi and the cone Pi. One 

o 

deduces the existence of an eigenvalue fio > of Sa, to which corresponds the eigenvector uq gPi, 
that is 

P^^+^uo = fiouo. (91) 

o 

Moreover, if ui gPi is any other eigenvector of Sa, related to a positive eigenvalue, then ui is 
equal to Uq up to a multiplicative positive constant. 
Step 2 (Existence). 

Denote by Spr(PQ,) the spectral radius of Ba, where Ba is considered an element of =Sf(L|^). It is 
well known (see for example Section III. 6. 2 in [23]) that 

Spr(P,)= hm ||Pril^a2v 
Since Uq G and Uq ^ 0, using ( fPT]) ) one has: 



II nmOo+3) II 
I RmOo+3)|| ^ 11^" ^" 11-^1/ _ „m 
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which triggers 

Spr(i?J > > 0. (92) 

On the other hand, let us denote by P^j the (reproducible) cone of positive functions in L\^. Due 
to the Weak Maximum Principle (Lemma 13. ip . 

B^iPl,) C Pi,. 

We also have that Ba is a compact operator from L\.j onto itself (due to the compact embedding 
of Hli in L|f, see |2l]). 

Use now the weak version of the Krein-Rutman theorem (see [22]) for the operator Ba, the Banach 
space L\j and the cone P|^. It turns out there exists an eigenvalue /iq > of Ba and a corresponding 
eigenvector uq G — {0}, i.e. 

BaUo = fj'OUo- 

Clearly tto also belongs to Hlj. Moreover, any other eigenvalue yU G C of is such that 

|/i| < /iQ. (93) 

Let us now denote Xq = 1 / flo — a; Xq is clearly an eigenvalue of L related to the same eigenf unction 
Uq. It then follows: 



Set (f = M in the above equation, and since / uq > { remark that uq G (fi) by the obvious 

Jq 

embedding Lp,{Q) C L^{Q)), we deduce that Aq = 0. Then the following expression: 



1p = — Mo 

■uodx 



n 
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gives a solution of eq fll7al) . 
Step 3 (Uniqueness). 

Assume ipi and ip2 are two solutions to the problem ( 117aP -( fT7cll . Then ipi and il)2 are non- negative 
eigenvectors of operator L corresponding to the eigenvalue 0. This implies 

Bai^m = -^m, m = 1, 2 

a 



which gives 



1\ 



io+3 



Sa'ipm = \ ~j '^rn, rU = 1,2. 



We then obtain ipm £ and by ( (!90|) ) we also have ipm GPi,m = 1,2. Now by the uniqueness 

property of the strong version of the Krein-Rutman theorem, there exists r > such that ipi (x) = 

r-?/'2(x), Vx G fi. Now since / ipi = / ip2 = b we obtain r = 1 which proves the uniqueness. 

Jn Jq 

Remark 4.1. From inequality ( (1931) ) one can deduce, proceeding in a classical manner, that 

Re{\) > 

for any other complex eigenvalue of operator L . The eigenvalue is then the principal eigenvalue 
of the operator L. Moreover, what is quite remarkable is the fact that is the principal eigenvalue 
of L for any function k G {W^'°°{fl))'^ . 



5 Final comments 

We have offered a proof to the fact that the FENE dumbbell configurational distribution function 
diffusion equation - see the corresponding boundary value problem described in eqs (I10ap - fll0d|) - 



has unique steady state solutions. In doing so, we relied on the Krein - Rutman theory of elhptic 
operators. 
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There are several motivations for this work. In pLj asymptotic solutions for the probability 
density diffusion equation - valid for slow flows - are presented, but no proof for the existence 
of such solutions is offered. While in this work we have proved the existence of solutions to the 
diffusion equation for slow and fast flows (that is irrespective of whether the velocity gradient is 
"small" or "large" ) , the questions related to the convergence of explicit expansion solutions given 
in P (and in what functional space it occurs) are still to be addressed to the fullest. Moreover, 
we expect our results to further the work in finding asymptotic solutions valid for "large" velocity 
gradients, i.e. for fast flows, for the FENE dumbbell model; a solution in such a case is known 
only for rigid dumbbells (see |25j). 

Now, the elastic (or rigid) dumbbell polymer chain models are certainly crude representations 
of the real chains. That set aside, they do capture several - but not all - important features of vis- 
coelastic flows (e.g. shear rate dependent viscosity, first normal stress difference). Moreover, they 
owe a certain popularity among polymer scientists (from experimentalists to applied mathemati- 
cians) to their relative simplicity. More realistic models use different chain representations. For 
instance, Doi and Edwards [26] developed the so-called tube model for melt systems, that makes 
use of the de Genne's reptation ideas of anisotropic chain diffusion. Schweitzer and co-workers 
[27] developed a mode-coupling model in which the intermolecular structural constraints upon the 
motion of single macromolecules are modeled as a many body caging effect. Ngai and Plazek 
developed their own coupling model [SSUSH], very successful in predicting the thermo-rheological 
complexity. This being said, bead-spring or bead-rod chain models still attract significant at- 
tention: see for example [301 ED E2l ESI [Ml [35]. For sure their full capabilities are still to be 
uncovered. 
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